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SUMMARY 

The purpose of the study covered by this report was to determine 
the feasibility of using an on- board digital computer and strap-down 
inertial reference unit for attitude control of an advanced OAO space- 
craft. An integrated, three axis control law was assumed which had 
been previously proven stable in the large. The principal advantage of 
this type of control system is the ability to complete three axis reori- 
entations over large angles. This system, although apparently complex, 
has the effect of simplifying the overall system. This is because a re- 
orientation is a simple extension of a hold or point operation, i.e. , mode 
switching may be simplified. 

The study indicates that a system of this type is feasible and offers 
many advantages over present systems. The time for a reorientation 
using this type of control is considerably shorter than the time required 
using more conventional methods. The computer and inertial reference 
unit usedinthe study had the characteristics of systems presently under 
development. Feasibility, therefore is within the present state of the 
art, requiring only continued development of these systems. 

The basic system is not limited to OAO but may be adapted for 
other three axis stabilized spacecraft. 
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A COMPUTER CONTROL SYSTEM FOR AN ADVANCED OAO 


INTRODUCTION 

In recent years, the microminiaturization of electronic components 
has made possible very small, general purpose (GP) digital computers. 

It is within the present state of the art to include these computers on- 
board spacecraft. One of the primary uses for such a computer is data 
processing. Another important use is in conjunction with the attitude 
control system of the spacecraft. The advent of the small computer has 
made it possible to consider control laws for the spacecraft attitude 
control system that are relatively complex. All the effects of dynamic 
coupling and large angle nonlinearities may be considered at the start 
of the design procedure. Small angle limits need not restrict the con- 
trol system design. Thus an integrated, three axis control system may 
be designed in place of three, small angle, single axis systems. 

Although a great deal of work remains in this area Mortensen (ref. 1) 
and Meyer (ref. Z) have proposed three axis control laws. 

Meyer defines the direction cosine matrix relating the spacecraft 
axes to a fixed inertial reference as the output of the attitude control 
system. This matrix is multiplied by a reference matrix, and the con- 
trol law is formulated from the product or error matrix. The control 
law he assumes is shown to be asymptotically stable in the large. The 
attitude matrix is assumed to be known at all times. It may be deter- 
mined by some type of inertial sensor (star trackers, gyros). The 
matrix formulation and multiplication, and the control law formulation 
require a computer. 

In an earlier work Mortens en proposed a control law under the 
assumption of a particular formulation for the spacecraft kinematics. 

He used the Cay ley- Rodrigue s parameters. These are also known as 
Euler -Rodrigues parameters (by Roberson) and the Gibbs vector com- 
ponents (ref. 3). Assuming the parameters defining the body orienta- 
tion are known, the control using these parameters he shows to be 
asymptotically stable in the large. The parameters are not easily 
measured but may be determined by solving the first order differential 
equations in the parameters and body rates. This method is similar to 
Meyer’s method except the attitude matrix is never formed explicitly. 
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The control law is formed directly from the kinematic parameters. The 
computation of the kinematic parameters and the control law formula- 
tion require a computer. 

Control laws of these types have three characteristics which should 
prove very valuable in future spacecraft: 

• Stability for any attitude. Thus the desired attitude of a space- 
craft is unrestricted and may be determined from considerations 
other than stability. 

• Three axis reorientation capability. Slewing all three axes 
simultaneously results in a rapid reorientation. 

• Simplicity. Since a three axis reorientation is simply an exten- 
sion of a hold or point operation, mode switching may be simpli- 
fied. 

This report is a partial result of a very short conceptual study for 
an advanced OAO spacecraft. The spacecraft considered was similar to 
the present OAO but with an extended telescope. The basic assumptions 
underlying the study with respect to the attitude control system were the 
presence of a GP, on-board computer and a very precise, strap-down 
inertial reference unit (1RU) 

One of the objectives of the study was the minimization of the num- 
ber of gimballed star trackers. The results of the study show that with 
the type of control law covered in this report the control system requires 
no gimballed star trackers. The trackers assumed are: 

• Boresighted star tracker along the roll axis parallel to the ex- 
periment telescope. 

• Fixed Canopus tracker along the spacecraft yaw axis for initial 
stabilization. 

• Two axis sun tracker, gimballed in pitch, for initial stabilization 
(See Figure 5). 

The primary purpose of the portion of the suidy covered by this re- 
port was to determine the feasibility of a three axis control system for 
pointing and reorienting the spacecraft. Assuming utilization of the 
computer and IRU, implementation of the type of control law proposed 
by Mortensen was considered. 
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It should be noted that although an advanced OAO is the specific 
spacecraft considered, many of the principles may be applied to other 
spacecraft that require a pointing and slew capability. 


CONTROL LAWS 

As in any spacecraft with a broad operating regime (initial stabili- 
zation through fine pointing) a number of different operating modes are 
required. This is true primarily because of the dynamic range limita- 
tions of the sensors and actuators. The control modes may be considered 
in three categories: 

• Initial stabilization. The initial three axis stabilization of the 
spacecraft requires the acquisition of two celestial references. 

The two references in this case are the sun and the star Canopus. 
Initial stabilization takes place in two stages. First, the sun line 
is acquired and two axis control about this line is maintained with 
the sun tracker. Acquisition of the sun line with a sun tracker and 
rate gyros (or IRU) requires a very simple control law. This con- 
trol law may be shown to be globally stable (ref. 3). No computer 
is required to implement it. In the second stage of the procedure 
the spacecraft is slewed around the sun line until the Canopus 
tracker detects and locks on Canopus. 

• Hold or Pointing. Operation in this type of mode, with the sun 
tracker, the experiment telescope, control system star tracker 
or gyros (IRU), is essentially small angle control. A linear 
position plus rate control [6 t kd), using either a lead network 
or gyro information, will give satisfactory performance. No 
computer is required. 

• Reorientation. This may be accomplished by three consecutive 
single axis slews. However, this is time consuming and com- 
plex. A large angle, three axis reorientation which is fast and 
simple may easily be accomplished with a computer, IRU and 
control law as previously described. 

The computer control system encompasses the hold or pointing 
mode as well as the reorientation mode. Therefore, once the initial 
acquisition is complete a single control law may be used for holding, 
pointing and slewing. This control law would use the IRU as the basic 
sensor during a reorientation and while holding (e.g. when the experi- 
ment is occulted), and the optical trackers when pointing. 
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This control system may also be used in the Canopus search phase 
of the initial stabilization. 

The remainder of the report is concerned primarily with this type 
of control system. 


AN INTEGRATED THREE AXIS CONTROL LAW 
Introduction 


In order to effect a three axis reorientation it is necessary to know 
the spacecraft attitude at all times. This requires knowledge of the or- 
thogonal transformation matrix relating the spacecraft body axes to a 
known, fixed inertial frame. Explicit determination of the entire matrix 
is not necessary as long as the three independent parameters of the 
kinematic representation are known. 

One method of paramaterizing this matrix is by the Cayley-Rodrigues 
parameters. See equations A-17 in Appendix A. In terms of these three 
independent parameters the matrix* is (ref. 3): 

1 + a 2 - fi 2 - y 2 2 (a /3 + y) 2 (ay - fS) 

L l 2 (a/3 - y) 1 _ a 2 + /3 2 - y 2 2 (/3y + a) (1) 

1 + a 2 +/3 2 + y 2 

2 (ay + /3) 2(j3y - a) 1 _ a 2 - B 2 + 7 2 


Continuous knowledge of a,/3, >, determines body attitude at all times 
with respect to the inertial frame. The parameters are not easily meas- 
ured directly but may be determined by continuously solving the differ- 
ential equation relating a,/ 3 , yto the body rates co * , , co z . 


1 + a 2 

a/3 - y 

ay + p 

a/3 + y 

1 +/3 2 

By - a 

ay - /3 

( 3y + a 

1 + y 2 


( 2 ) 


^Capital letters with double overbars as square matrices; with single overbars are column matrices. 
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Using these parameters Mortensen has proposed the following con- 
trol law (ref. 1): 


u x - k 2 i w x + k i (1 + °- 2 + @ 2 + y 2 ) a 

u y = k 22 0) y + kj (1 + a 2 + /3 2 + y 2 )fi 

u z ~ k 23 a) z + hj(l + a 2 + /3 2 + y 2 )y 

where u x , u y , are control torques. 


(3) 


He shows the system with this control law to be asymptotically stable 
in the large with respect to the origin (a=/3= : y=0). The dynamic and con- 
trol equations are developed in Appendix A. This control law is analogous 
to a conventional law in that the torque is a function (although non-linear) 
of rate and position. For small angles (i.e. the spacecraft body axes near 
the desired attitude) the following approximations hold (ref. 3): 
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a ^ 


2 


(4) 


y 


where d x , d r 6 3 are Euler angles. 


Therefore, for small angles the control reduces to a simple linear 
rate plus position control. 


u 
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y 




k 22 M y 


+ k 3^ 2 


U z ~ k 23 aJ r +k 3^3 


(5) 


A control law of this type, globally stable for any attitude, encom- 
passes pointing and slewing without requiring mode switching. The 
inertial reference is defined as a = j3-y = 0. When a new reference is 
desired e.g. pointing at a different star, the transformation relating the 
present position to the new position is determined. The corresponding 


5 



parameters a,/3,y are determined and sent up to the spacecraft. The 
spacecraft, utilizing this control law, acts to null a,/3,y and thereby 
moves to the new attitude. This is a reorientation or three axis slew. 

Using the strap down IRU and GP computer, a ,8$y may be continu- 
ously updated to allow utilization of this control law. See Figure 1. 


Attitude Reference Algorithm 


The sensor to be used for operational control of the spacecraft dur- 
ing a reorientation is an IRU composed of three inertial gyros. This is 
a strap-down system with the gyros operating in a pulse rebalance loop. 
They are continuously nulled by a series of pulses the size and rate of 
which provide position and rate information. References four and five 
give a good description of this type of gyro. The information from the 
IRU must be used in conjunction with the computer to update the param- 
eters a,/3,y throughout a reorientation. This may be done by solving 
equation 2. 

The body rates are not known accurately enough for equation 2 to 
be solved in its present form. A similar problem was encountered in 
reference five with respect to the updating of an attitude reference 
matrix. An integration algorithm based on an incremental angle was 
developed. A similar procedure may be used here. 

A Taylor series expansion may be used to update a>/3$y . The value 
at time t+h may be determined from the value at time t and the informa- 
tion received from the gyro during the sampling interval h. Only the 
parameter a(t) will be considered here although similar expressions 
exist for /3(t) and y(t). 


a(t + h) = a(t) + a(t)h + i a(t) h 2 + \ a(t) h 3 + . . . (6) 

2 6 

Using only a first order expansion and substituting for a(t) from equa- 
tion 2: 

(1 + a 2 ( t)) A (a(t) yfl(t) -?(0) A y (a(t) 7(t) +/3(t)) A 
a(t + h) = a(t) + i + - M?) 

where A. = co. (t)h; i = x, y, z. 
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Figure 1-Control System Block Diagram 













The A. terms represent the angular output of the i th gyro during the 
interval h, i.e. a number of pulses of weight q (q is the gyro quantiza- 
tion level, 2.4 arc sec in the IRU slew mode). The updating therefore 
may be done on an incremental angle basis. If the second order terms 
of the Taylor series are kept, terms on the order of h 2 appear. Ap- 
pendix B lists the second order expansions for all three parameters. 

Updating the parameters in this manner allows the control law to 
be updated at each sampling point also. The body rates used in the con- 
trol law are simply Ai/h. 

A study was conducted to determine the feasibility of this type of 
algorithm. A digital computer simulation of the system was used to help 
evaluate the errors. The attitude error was defined as the difference 
between the true spacecraft attitude and the attitude computed by means 
of the Taylor expansion. The spacecraft dynamics and true attitude were 
computed by means of a four point Runge-Kutta routine. The integration 
step was kept much smaller than the sampling interval h. 

The results of the study are as follows: 

• The error in the update algorithm does increase with an increase 
in the sampling interval h. This is not a serious constraint be- 
cause even relatively long sampling intervals gave good results. 
Using a second order Taylor expansion and a 5 second sampling 
interval the errors at the end of a 60° slew were on the order of 
14 arc seconds. 

® The error in the update algorithm does increase with the mag- 
nitude of the slew angle for the longer sampling intervals (2 sec- 
onds and up). For a sampling interval of 100 ms the error is 
apparently independent of the magnitude of the slew angle even 
for very long slews (165°). It remains on the order of 2.4 arc 
sec which is the quantization level of the gyros. For a sampling 
interval of one sec the error for a 165° slew is only about 5.3 arc 
sac. This was using second order Taylor expansion. 

• For moderate sampling intervals and slew angles (up to 2 sec- 
onds and 90°) the attitude error is primarily a function of the 
gyro quantization level and is on this order of magnitude. This 
also holds true for shorter slews (30°) with longer sampling 
intervals (5 seconds). 
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• The computer round-off error may be kept negligible by using a 
36 bit, double precision word. The computer considered in the 
study (see the following section) has the ability to operate single 
(18 bits) or double precision. 

The estimated times for the update computations are shown in 
Table 1. 


Table 1 


Estimated Time for the Update Computation 



Single Precision 

Double Precision 

1 st order Taylor 

2 ms 

7 ms 

2nd order Taylor 

4 ms 

14 ms 


These times are based on the computer described in the following sec- 
tion. These computation times are very compatible with the sampling 
intervals investigated. 

These results indicate this type of algorithm for updating the space- 
craft attitude parameters is feasible. 


Computer 


A candidate for the computer required for attitude updating and con- 
trol law formulation is the "Units” On-Board Processor being developed 
by the Space Electronics Branch, Information Processing Division, GSFC 
(ref. 6). This computer is being developed primarily as a data processor 
but appears to have the capacity and capability for handling the control 
equations. Some of the features of this computer are: 

• Post launch reprogramming capability. The control laws may 
be changed in flight if desired. 

• Plated wire memory unit with 8192 word (18 bits/word) capacity. 
This may be easily expanded to about 65,000 words. 

• Low Power. Total power required is about 10.7 watts with peak 
memory activity and less than one watt at stand-by. 

• Very fast, Duty cycle of less than 0.Z5 at a Z/xs write rate. 
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• Very small. Memory size is about 3X6X10 in. (For 8192 words). 

• Maximum use of monolithic integrated circuits, 80% utilization. 


Inertial Reference Unit 


The inertial reference unit (IRU) being developed by MIT for OAO 
will apparently meet the attitude reference requirements. Three opera- 
tional modes are available: hold, slew and I-stab. The primary concern 
of this report, three axis reorientation, requires use of the slew mode. 
For this study the processing being developed for the early models of the 
IRU was not considered. The only concern was the basic performance of 
the gyros. 

Following are some of the important characteristics of the IRU in 
its slew mode: 

• Resolution, 2.4 arc sec/pulse 

• Maximum rate, 480 arc sec/sec 

• Pulse rate, 200 pps 

• Compensated drift, 20 arc sec/orbit (12 arc sec/hr) 

• Torquing error less than 0.01% 


Normal Reorientation 

The system dynamic, kinematic and control equations were simu- 
lated on a digital computer. The primary objective of the simulation 
was to determine performance with respect to a slew or reorientation. 
The equations reduce to approximately linear equations for small angles 
(hold or pointing). Thus it was not necessary to consider the hold or 
pointing modes of operation. 

The angle <fc (see Appendix A) is a good scaler representation of 
error, i.e., difference between the actual body orientation and the de- 
sired orientation (ref. 2). This is true because a , /3 , y are all zero 
if and only if <£> is zero. Therefore was chosen as the primary per- 
formance index. For a single axis slew 4> is equivalent to the Euler 
angle around the slew axis. 
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The only requirements for stability are that the position gain (k p ) 
and rate gains (k x , k y , k 2 ) be positive and non-zero. There are no other 
requirements with respect to magnitude or linearity. Therefore once 
the stability criteria are met the designer may set the gains to meet 
other criteria such as speed of response, damping, etc. 

Response . Probably the most important criterion for a normal re- 
orientation is the time for a slew. The shorter the reorientation time, 
the longer the experimenting time. Another criterion, maintenance of 
a stable slew axis, is discussed in a following section. Using the time 
criterion several computer runs were made. 

The time for a slew was defined as the time required from the read- 
in of initial conditions to the on-board computer until the norm reached 
its minimum value. The norm; defined as 


NORM = + (or 2 + a; 2 + cc 2 + 0 2 ) 1/2 

x y z 


( 8 ) 


was usually assigned a minimum of 10“ 4 . Since the norm is always 
positive, the rates were required to be below 10 4 rad/sec and <£> was 
required to be below 10” 4 rad (0.34 arc-min) before a slew was con- 
sidered complete. Therefore slew or reorientation time as defined here 
includes settling time. 


Since a slew is basically a momentum exchange the gains were set 
to result in a maximum wheel speed during the slew. The rate gains 
were set to provide good damping at the end of the slew. Both type A 
and B wheels were used. See Figure 2. The characteristics were as 
follows : 


Stall Torque 
Max. Momentum 
Time Constant T 

m 


0.27 Nm (0.20 ft-lb) 

13.6 Nms (10.0 ft-lb-s) 

50 s 


To compare a three axis reorientation with the more conventional 
method using three single axis slews, a number of runs were made 
with the slew axis eigenvector aligned with a control axis. This is 
equivalent to a conventional single axis slew. 


Figure 3 shows the response for four single axis and two multi- 
axis slews (eigenvector 0.5, 0.5, 0.707). The maximum slew rate for 
the single axis slew is about 0.25 X 10~ 2 rad/sec. This is determined 
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by the 13,6 Nms maximum wheel momentum and the 5420 kg-m 2 body 
inertia (the same for all axes). The maximum slew rate for the multi- 
axis slew occurs when all three wheels are saturated. It is about 0.43 X 
10“ 2 rad/sec. Thus all three wheels are being used to reduce 0. The 
single and multi-axis slews in Figure 3 are not directly comparable be- 
cause the final orientations are different. A comparison will be made 
later in the report. 

Figure 4 is a plot of total slew time for various length slews. The 
curves all converge to about 160 sec for a 2° slew. This is apparently 
due to the end effects, i.e., acceleration at the start and settling at the 
completion of the slew. The slope of the linear portion of the single 
axis curves is (0.25 X 10~ 2 ) -1 sec/rad. This is determined by the maxi- 
mum wheel momentum. All the points checked for multi-axis slews were 
below the time required for equivalent length single axis slews. 

To compare the time required for reorientation by means of a multi- 
axis slew with the time required for the same reorientation via three 
single axis slews, three points were checked both ways. See points A, 

B, C on Figure 4. Table 2 summarizes the results. 


Table 2 

Reorientation Times 



Euler Angles 
(rad; deg) 

Eigenvector Comp. 

4> 

(rad; 

deg) 

T, 

(sec) 

T 

1 3 

(sec) 

T / T 

13 

e i 

* 2 

*3 

c 

X 

C 

y 

C 

Z 

A 

0.1745 
10 . 

0.2745 

15.7 

0.1745 

10 

0.5187 

0.6797 

0.5187 

0.3797 

21.8 

700 

260 

2.69 

B 

0.523 

30 

0.523 

30 

0.523 

30 

0.6546 

0.3782 

— 

0.6546 

0.9722 

55.7 

1080 

420 

2.57 

C 

1.045 

60 

j 

1.045 

60 

1.045 

60 

0.6946 

0.1869 

0.6946 

2.0327 

116.5 

1680 

790 

2.13 


In the table T x is the total time for three single axis slews. This 
time is determined by using the Euler angles and single axis slew curves 
of Figure 4. The time to complete the same reorientation using a multi- 
axis slew is T 3 . The ratio Tj / T 3 shows the multi-axis reorientation 
to be over twice as fast as the conventional slew. This should be ex- 
pected for at least two reasons: 
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• The end effects which are independent of length of slew occur 
three times for three slews and only once for a single slew. 

• All the wheels contribute simultaneously during a multi-axis 
slew causing the total rate which the error is being reduced to 
be higher than for a single axis slew. 

The slew time using a type B motor was slightly longer in each case 
due to the reduced torque available for acceleration and deceleration. 

The slopes of the linear portions are the same since the maximum wheel 
speeds are the same. 

Accuracy - In a hold or point mode the accuracy is primarily a 
function of the sensors. When the spacecraft is holding on the experi- 
ment telescope or roll axis star tracker there is no problem. When 
holding on the inertial reference unit (IRU) the primary error will be 
due to the gyro drift. This is expected to be about 12 arc sec/hr. The 
desired accuracy determines the update frequency when the spacecraft 
is holding on the IRU. 

There are three sources of error present during a reorientation. 
They are the attitude updating, gyro drift and the torquing inaccuracy 
of the gyros. 

For a sampling interval of 100 ms the attitude update error may be 
kept on the order of the gyro quantization level which is 2.4 arc sec. 

The reorientation error due to gyro drift will be small because 
most slews will take less than 20 minutes to complete. The compen- 
sated drift during this time is only 3 to 4 arc sec. 

There is an error accumulation due to the torque rebalance loop 
when the gyro is being torqued. The torque rebalance pulses will con- 
tain some error in size and shape. This will cause an attitude error 
build-up over a long slew. The torquing inaccuracy is estimated to be 
0.01%. This represents an error of 0.0165° (one arc min) over a 165 
slew. Therefore, this is the major attitude error in a large angle 
reorientation. 

The accuracy of the system need not be as high as that required for 
fine pointing. The objective of the reorientation is the acquisition of a 
new experiment star. Therefore, the primary factor in determining the 
required slew accuracy is the field of view of the experiment telescope 
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or the roll star tracker. If the total reorientation error is well within 
the field of view of the sensor it is feasible to reorient the spacecraft 
by the method described. With an experiment that has a field of view 
±4 arc min (e.g. OAO-Pr inceton Experiment Package) it does appear 
feasible to move from fine pointing to fine pointing via a single, large 
angle, three axis reorientation. 

If the experiment field of view is much smaller, it is possible to 
move to a 6th magnitude or brighter star and lock on this star with the 
roll star tracker. The IRU is reset and the spacecraft is then moved to 
the experiment star. This is practical because the average distance 
from any random point in the sky to the nearest 6th magnitude or brighter 
star is only about 1.5° (ref. 7). See the section on updating. The error 
for this magnitude slew will be on the order of 2.5 arc sec. The time 
for the short slew will be about two minutes. See Figure 4. 

Disturbance — The effect of environmental disturbance torques 
during a slew is to increase the total momentum of the system. The 
reorientation trajectory with disturbances will in general be different 
from the undisturbed trajectory. However, the final reference as stored 
in the computer and IRU is not affected; the system will null at the same 
reference as the undisturbed system. The time for the reorientation 
will depend on the magnitude and direction of the disturbances and may 
be shorter or longer than that required for an undisturbed reorientation. 
Although an extensive torque disturbance analysis was not conducted 
some computer runs were made with disturbance torques. No gross 
changes in system performance were noted. The additional time re- 
quired for one 165° slew was on the order of one minute with disturbances 
on the order of 0.0001 to 0.0005 Nm (1000 to 5000 dyne-cm). 


Canopus Search 

As stated previously the time for a slew is usually an important 
criterion. It is, however, not the only one. For a Canopus search 
maintenance of a constant slew is more important. 

In the sun, a Canopus search by slewing around the sun line is 
easily accomplished. Position control with respect to the sun line is 
maintained by the sun tracker. The roll and yaw wheels are biased as 
a function of the angle between the sun line and roll axis. See Figure 5. 
Due to the time of year, orbit inclination, etc. it may not be possible 
or advantageous to do a Canopus search in the sun. It then becomes 
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necessary to slew around an occulted sun line. In order to insure that 
the fixed Canopus tracker does acquire Canopus, the slew must remain 
fixed within about 0.2° during the slew. This is equivalent to a three 
axis reorientation as described previously but with the spacecraft slewing 
around a fixed eigenvector. This may be accomplished by commanding 
a reorientation around this eigenvector and limiting the motor torques 
and maximum momenta by the ratio of the eigenvector components. The 
body rates are proportional to the eigenvector components and the total 
velocity vector remains collinear with the slew axis or eigenvector 
throughout the entire slew. The Canopus tracker trajectory is therefore 
well defined. 

In the normal reorientation the location of the eigenvector is un- 
constrained. For a significant portion of the slew all three wheels are 
saturated and the body rate vector is not collinear with the eigenvector. 
Because of the fixed relationship of the body rates a constrained rate 
slew will be slower than a normal reorientation. However, if at least 
one wheel is allowed to saturate, the slew rate will be higher than that 
for a single axis slew. 

A computer simulation confirmed the stability of the eigenvector 
position for this type of reorientation. 

This type of Canopus search is essentially a position reorientation 
around the occulted sun line. A large slew angle (up to 180 ) is com- 
manded around the desired slew axis. When a Canopus presence signal 
is sensed by the Canopus tracker the IRU is reset and the mode switched 
to IRU-computer control for settling. After settling, roll and pitch are 
available from the Canopus tracker and yaw is held by the yaw gyro 
relative to the sun line. Once in the sun the IRU may be switched from 
its I-stab mode, reset and retrimmed. The gyro quantization level in 
the I-stab mode is 38 arc sec. With a 1° FOV Canopus tracker this 
level should be adequate to insure Canopus comes within the tracker FOV. 

Because of the limited maximum slew angle it may be necessary to 
command a second or at most third slew. However, an attitude reference 
with respect to the non slew axes is maintained throughout the entire slew. 


Updating 

The updating of the IRU requires some type of celestial sensor 
(star and/or sun trackers) to relate the spacecraft axes to the celestial 
reference. Updating consists of two operations, a reference angle reset 
and a drift biasing or retrim. 
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With the boresighted star tracker (BST) along the roll axis, resetting 
the pitch and yaw gyros is relatively simple. These gyros may be reset 
any time the spacecraft is holding on the experiment or roll BST, 

The updating of the roll axis gyro is slightly more complex because 
of the tracker locations. A side looking tracker may be used, such as 
the gimballed sun tracker or an electronically gimballed star tracker. 
Both methods have advantages and disadvantages with respect to field 
of view, gimbal reliability, location of stars, etc. These will not be 
discussed further here. 

The relative speed and accuracy of a three axis reorientation makes 
it possible to reset the roll gyro by another method using only the roll 
BST. The equations are developed in reference 8. The procedure is as 
follows : 

1. With the spacecraft holding on a known star with experiment or 
roll BST, the pitch and yaw gyros are reset. The direction of the roll 
axis in inertial space is now known but some small unknown roll angle 
exists. 


2. Reorient the spacecraft so the roll axis is pointing at a second 
known star. A two axis pitch/yaw error will exist in the roll BST due 
to the initial roll error. The initial roll error must not be so large 
that the second star is out of the field of view of the roll BST at the end 
of the slew. 

3, With the known transformation matrices and errors read from 
the roll BST compute the roll error and reset the roll gyro. The compu- 
tation may be done on board or on the ground. 

The primary advantage of this method is its independence of a side 
looking tracker for a roll reset. 

Because an extra slew is required with this method, the separation 
of stars that may be used by the roll BST is important. A study was 
carried out to determine the distribution of 6th magnitude and brighter 
stars. Starting at a random point in the sky the distance to the nearest 
unambigous 6th magnitude or brighter star was determined. Figure 6 
is a plot of the data for 1000 random points. An unambiguous star is one 
that has no other stars of the same magnitude or brighter within 14 
arc -min of it. A BST with a 10 arc min square field of view will see 
only one of these stars. Of 5300 stars (6th magnitude or brighter) about 
5000 are unambiguous. The average distance that must be moved is about 
1.5°. 
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DISTANCE TO NEAREST STAR (0E6.) 


The time to complete a slew this short is only about 2 minutes. 
Therefore, the time for a reset on all three gyros would be short. Even 
if a reorientation was necessary to acquire a star for the pitch and yaw 
gyro reset, the total time required for a reset would generally be less 
than 10 minutes. 

A retrim or drift biasing requires a hold on a tracker for the time 
required to check gyro drift. Retrimming therefore does require a side 
looking tracker. * Since retrimming would not be required as often as 
resetting, this method of resetting, independent of the side looking tracker 
offers systems flexibility. 


DISCUSSION 

This study was not a general study of various control laws and/or 
kinematic representations. The choice of Mortensen's control law with 
the Euler-Rodrigues parameters appeared to be natural for this problem. 
Further study may show it is advantageous to use Meyer's method where 
the entire transformation matrix is updated. Basically there is little 
difference between the two approaches. It may be shown the control laws 
are very similar. 

These types of control laws may be used with other types of sensors. 
In particular the attitude matrix defined by Meyer as the position output 
of the system is independent of the sensors used to describe its elements. 
Star trackers may be used as well as gyros. However it appears that a 
strap-down inertial reference unit as considered here offers the most 
advantages for large angle reorientations. The use of gimballed star 
trackers has two disadvantages. Due to the limited gimbal rotation, 
large angle reorientations cannot be completed without switching trackers 
during the slew. This introduces additional constraints on the reorienta- 
tion maneuver. Secondly, the attitude matrix is more complex because of 
the presence of the trigonometric functions resulting from the tracker 
gimbals. 

Because of the gyro drift it appears necessary that the spacecraft 
hold its position during the experimenting time with optical trackers, 
either the experiment telescope or the roll axis boresighted star tracker. 
Holding the spacecraft position with an optical tracker requires a 


*Roll retrim is theoretically possible with the roll BST by commanding equal and opposite pitch 
slews separated by the measuring time interval. Using the roll update technique the roll drift 
during the interval can be calculated. 
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relatively simple control system and does not require a computer. 
Therefore, during this time the computer may be freed for processing 
the experiment data* During a reorientation when no experimenting 
is being carried out, the computer may be used completely for the 
control system. Although the computer controlled system may easily 
handle both the hold and reorientation mode it may be more efficient 
to use it in this manner. 

The in-flight programming capability of the computer may be used 
to advantage* The control gains may be set initially to carry out the 
constant eigenvector slew that is required for a Canopus search. Once 
this is complete the gains may be changed to allow a more rapid slew 
for normal reorientations. 


CONCLUSIONS 

1. A computer control system that performs as well during a slew 
as during pointing is feasible for an advanced OAO. Utilizing an on- 
board digital computer and strap-down inertial reference unit the 
system may complete large angle, three axis reorientations with the 
same control law used for pointing. The computer speed is such that a 
second order Taylor series may be used to update the spacecraft attitude 
during a reorientation. The errors due to the attitude updating may be 
kept on the order of the gyro quantization level. The primary error in 
a large angle reorientation is due to the gyro torquing inaccuracy. This 
error, however, is not prohibitive. 

Z. A general purpose computer now under development appears to 
meet the computational requirements. This is the M Units n On-Board 
Digital Processor being developed by the Space Electronics Branch, GSFC. 

3. An inertial reference unit being developed for GSFC by MIT 
appears capable of meeting the sensor requirements. 

4. A three axis slew is superior to three single axis slews when 
time is the reorientation criterion. A three axis slew was shown to 
be over twice as fast as three single axis slews for the reorientations 
that were checked. It is also simpler, requiring less commands at the 
start of the slew. 

5. A Canopus search on the dark side of the orbit may be carried 
out with the same type of control law. Using the proper gain settings 
a slew may be maintained about a constant, arbitrary axis established 
by the sun tracker. 
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6. The inertial reference unit reference may be reset with a single 
boresighted star tracker. This may be accomplished by successively 
pointing at two separate stars with the roll axis boresighted star tracker. 
The complete resetting procedure would take less than 10 minutes. 

7. It is possible to use a conventional linear rate and position 
control when pointing and employ the computer primarily for reorienta- 
tion. The computer could then be used for data processing when the 
spacecraft is in a pointing mode. 


REFERENCES 

1. Mortensen, R. E., On Systems for Automatic Control of the Rotation 
of a Rigid Body, Electronics Research Laboratory, University of 
California, Berkeley, November 27, 1963. 

2. Meyer, G., On the Use of Euler’s Theorem on Rotations for the 
Synthesis of Attitude Control Systems, Ames Research Center, 

Moffett Field, California, NASA TN D-3643, September 1966. 

3. Sabroff, A., et. al., Investigation of the Acquisition Problem in Satel- 
lite Attitude Control, Technical Report, AFFDC-TR-65- 1 1 5, Air 
Force Flight Dynamics Laboratory, Air Force Systems Command, 
Wright Patterson Air Force Base, Ohio, December 1965. 

4. Otten, D. D., A Look into Strap-down Guidance Design-I, Control 
Engineering, October 1966. 

5. Otten, D. D., Body Fixed, Three Axis Reference System Study - 
Phase I FinafReport , TRW Systems, TRW No. 4499-6007-R0000, 

May 2, 1966. 

6. Taylor, T., et. al., A General Purpose On Board Processor for 
Scientific Spacecraft, GSFC X-562-67-202, July, 1967. 

7. Fang, A. C., Determination of the Proximity of Acceptable Reference 
Stars for Arbitrary Spacecraft Orientations, Stabilization and Control 
Branch Report No. 161, GSFC, June 1, 1967. 

8. Hrastar, J. A., A Digital Control System for an Advanced OAO, 
Stabilization and Control Branch Report No. 163, GSFC, June 9, 1967. 


24 



APPENDIX A 


The spacecraft dynamic and control equations are developed in this 
appendix. Although Mortensen did not use reaction wheels for control, 
the method for proving stability closely parallels his development (ref. 1). 

Assuming the control torque on the spacecraft is developed by 
reaction wheels the dynamic equations are (ref. 3). 


H + (I - I ) co co +(ojH - co H )=M - H 

x '• z y'yz ^y wz z wy ' x wx 

H +(I - I )co co +(wH -a) H ) = M - H 

y v x z ' x z v z wx x wz' y wy 

H -f (I - I ) co co + (co H - co W )-M - H 

z V y x ' . x y v x wy y wx ' z wz 


where 

H x , H y , H z — spacecraft momentum about x, y, z (principal 
axes) less wheel momentum. (Also H i = I £ ; 
i = x , y, z). 


I , I , I - spacecraft inertia about x, y, z. 

x y z 


H wx , , H wz - wheel momentum about x, y, z. 


M x , M y , M z — external torques about x , y , z . 


Assume M x = M y = M z = 0. During a slew the external torques are 
small compared to the control torques and the assumption should be 
valid. Putting the equation into matrix form: 



H 

X 


0 

-Hrz/Iy 

Hry/I. 


H x 


H wx 


H 

y 


Hr,/ 1 . 

0 



H 

y 


Ky 


k 


_->V I x 

Krx/I y 

0 


H 

z 


k. 

where H . = 

Ti 

• H. 

i 

+ H . ; 

W1 

i = x, y, 

z. 





(A-2) 
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Combining equation A- 2 with the kinematic equation 2: 



Assuming a control law of the type shown i: 
wheel with torque -speed curve as shown in 
equations are: 



l equations 3 and a reaction 
Figure 2, Type A, the control 
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The total system is ninth order requiring nine state variables for a 
complete description of the trajectory in state space. The state variables 
are the three body momenta, three parameters specifying body position 
and the three total momenta. The equilibrium position of the spacecraft 
may be shown to be independent of the system total momenta. Thus 
there is a three dimensional subspace every point of which is an 
equilibrium point for the spacecraft. Therefore, for spacecraft stability 
analysis the six dimensional vector representing the body position and 
' rate may be used in place of the nine dimensional system vector. The 
total momentum in this formulation may be considered a time varying 
parameter. 


H x 


-k 

X 

-Hr/Jy 

«t/i. 

M'f) 

H 

y 


h Tz /i x 

-k 

y 

-Ktx/I, 

-k( 

'a/3\ 

2 ) 



-Hry/ 1 . 

H Xx/Iy 

-k 

z 

-k 

p 

(¥) 

a 


1+a 2 

21, 

a f3 

2I y 

a y 

2I z 


0 

b 


a /3 

2 I 

X 

1 +/3 2 

21 

y 

fiy 

2 I 

z 

-H./2I, 

y 


a y 

f3y 

1 +y 2 

H /2I 



2 I 

X 

2 I 

y 

2 I 

Z 


y j 


This is in the form 

X = F (X) X 

where X - [H x , H y , H z , a, /3, y] T 

F (X) — nonlinear system matrix. 





H /2I 

z' z 

0 

-H/2I 

x' X 



-H y /2I y 

V2I X 


0 


(A-6) 
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Define the matrices K and R as follows 


K = dia g [ -k x , -k y , -k z , 0,0,0] 
R = F - K 


Choose as a candidate V function the quadratic form 


V = X T G X 


where 


G = diag. 

Since 0 is a constant it follows that 


* * 1 k k k 

T’ T’ T~' p’ p* p 
L x y * 


v = x t gx + x t gx 


Substituting equations A- 7 into A- 6 and then into A-9 it (A 

v--x T r t gx + x t k t gx + x t grx + x t gkx 

Note that G R is skew symmetric and that G = G T , therefore 

G R = - (G R) T = - R T G t = - R T G 
This reduces equation A- 10 to 

V = X T K T G X + X T G K X 
The following also holds 

K T G = K T G T = (G K) T = G f 


(A-7) 

(A-8) 

(A-9) 

9) becomes 
(A- 10) 

(A - 1 1) 

(A- 12) 

(A-13) 


28 



Therefore 


Thus if 


V = 2 X T G K X 



k > 0 , k >0, k >0 

x y z 


then V is negative semidefinite. Since 


V = 





a 2 + k p /3 2 


+ k P 


7 2 


(A- 14) 


(A- 15) 


(A- 16) 


is positive definite (k p > 0) the origin is stable. It may be shown that 

the origin is not only stable but asymptotically stable, i.e., V - 0 as 

t - oo . 

Since following relations are valid (ref. 3) 

a = C x tan 0/2 

ft = C y tan 0/2 (A- 1^) 

y = C tan 0/2 


where 


C , C , C - the components of the +1 eigenvector of the orthogonal 

x y z 

transformation parameterized by a , ft, -y 

4> — the angle around the +1 eigenvector that the body frame 
is rotated from the inertial frame. 


equation A- 16 may be written 




h 2 

+ — + k p tan 2 4>/2 


(A- 1 8) 
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Therefore, for V to be bounded it is necessary that the body angular 
momentum be bounded and 0 must be less than rr . The total system 
momentum must be within the total capacity of the wheels. It will be 
stored in the wheels at the completion of the reorientation. The latter 
constraint (0 < n) is not a practical restriction on the system since it 
is easily satisfied. With the boundedness of V established it may be 
concluded that V is a valid Lyapunov function and the stability con- 
clusions are valid. 

Asymptotic stability in the large has not yet been proven in general 
when a type B reaction wheel is used. However two important special 
cases may be shown to be asymptotically stable. One is the case in 
which the total system momentum is zero. The other is the linearized 
system near the origin. 

When a type B reaction wheel is used the additional back emf terms 
- H wi / T m i i = x , y , z are added to the respective control equations 
(A-4). When the total system momentum is zero H. = - H wi • Thus the 
system is reduced to sixth order and it is easily proven stable in a 
manner analagous to that used in the preceding proof. 

When the system equations are linearized near the origin the result 
is three, uncoupled, linear, second order systems. These are easily 
proven stable by conventional methods. Any stored momentum repre- 
sents a shift to a new origin. 
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APPENDIX B 


The equations for updating the parameters a, /3, y using a second 
order Taylor series are the following: 

a „ + l = a n + J £(1 + a n) A x + ( a rA " ?„> A y 
+ (a n7n+ /3n) A J > A x 

+ < a n " ?„> A y + >n < a n ^ A z 

+ (2 a 2 /S +/3 -a 7 ) A A 

v n 'll “n n x y 

+ (2a 2 y +y + a ft ) A A 

v n / n / n n 7 x z 

+ ( 2 a n 4 +^n " A y A z^ 

4, + i = 4, + J [ K/? n + 7 n ) A x + (1 + /3 n 2 )A y 

+ >n - a n) A z^ + \ K ( a „ A x 

+ /3 n (i + /3 2 ) A 2 + y n 0 n y n - a n ) A 2 
+ (2a /3 2 + a +/3 y ) A A 

v n 'n n 'll 7 n 7 x y 

+ (2 a B y - a 2 + y 2 ) A A 

^ n 'll 7 n n 7 n 7 x z 

+ (2 /? 2 y +7 -a yS ) A A] 

v n 7 n 7 n n 'll 7 y z J 
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>n + l = + \ [ ( a n " 4.) A x + (4. + a n ) A y 

+ (1 + y 2 n ) A z ] + \ [a n (a n 7 n - /3J A* 

+ 4. (4, >n + a n) A y + 0 + ?n> A z 

+ ( 2 a n 4, + a n " ^n> A x A y 

+ ( 2 a 7 2 + a - /3 7 ) A A 

v n / n n 'll ' n ' x y 

+ (2 /? 7 2 + /3 + a 7 ) A A ] 

T v ' n 7 n T ' n n / n / y z 
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